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Abstract. The main theoretical aspects of gravitoelectromagnetism ("GEM") are 
"q I presented. Two basic approaches to this subject are described and the role of the 

gravitational Larmor theorem is emphasized. Some of the consequences of GEM are 
5— i ■ briefly mentioned. 

: 1.1 Introduction 

H , 

The analogy between gravitation and electromagnetism has a long history. The 
similarity between Newton's law of gravitation and Coulomb's law of electricity 
naturally led to a gravitoelectric description of Newtonian gravitation. Moreover, 
on the basis of advances in electrodynamics in the second half of the nineteenth 
century, Holzmiiller pQ and Tisserand [2] postulated that the gravitational force 
exerted by the Sun on the planets of the solar system had an additional "magnetic" 
component. This extra force led to the precession of the planetary orbits; therefore, 
it could be adjusted in order to account for the excess perihelion precession of 
Mercury. Decades later, however, Einstein's general relativity provided a beautiful 
explanation of the excess motion of Mercury's perihelion in terms of a relativistic 
gravitoelectric correction to the Newtonian gravitational potential of the Sun [3]. 
Furthermore, general relativity, which is a field theory of gravitation, contains a 
gravitomagnetic field due to mass current [4]. Indeed, to bring together Newtonian 
gravitation and Lorentz invariance in a consistent field-theoretic framework, the 
introduction of a gravitomagnetic field is unavoidable. 
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According to general relativity, the proper rotation of the Sun produces a gravit- 
omagnetic field and the influence of this field on planetary orbits was first considered 
by de Sitter [5] and later in a more general form by Lense and Thirring [3]. The 
gravitomagnetic contribution to the excess motion of Mercury's perihelion turns out 
to be much smaller and in the opposite sense compared to the main gravitoelectric 
motion; in fact, it turns out that the Lense-Thirring precession of planetary orbits 
is too small to be measurable at present. On the other hand, evidence for the grav- 
itomagnetic field of the Earth has been offered by Ciufolini by studying the motion 
of laser-ranged satellites LAGEOS and LAGEOS II [6]. The precise measurement 
of this field via superconducting gyroscopes in a drag-free satellite in polar orbit 
about the Earth is one of the aims of NASA's GP-B [7]. 

Within the framework of general relativity, gravitoelectromagnetism ("GEM") 
has been discussed by a number of authors [21 IS]; a more extensive list of references is 
provided in [9]. The purpose of this review is to present the two principal approaches 
to GEM and briefly describe some of their consequences. 



1.2 Linear Perturbation Approach to GEM 

We are interested in the general linear solution of the gravitational field equations [3]. 
It is assumed that a global background inertial frame with coordinates x^ 1 = (ct, x) 
and Minkowski metric tj^ is perturbed due to the presence of gravitating sources 
such that g^y = r]^ u + h^ u {x). It proves useful to define the trace-reversed amplitude 
h^u = h^ u — ^rj^h, where h = rj^h^u is the trace of h^ u . To linear order in the 
perturbation, Einstein's field equations 

1 8vrG 
R»v - -jQ^R — (1-1) 

take the form 

- 16-7rG m 

□v = — ( L2 ) 

after imposing the transverse gauge condition h jV = 0. The general solution of 
(jl.2p is a superposition of a particular solution together with the general solution of 
the wave equation; however, we are only interested in the special retarded solution 
of (|1.2p given by 

v=f/ T "V x ;r' l,x ' ) dV - m 

Let us define the matter density p and matter current j = pv via T 00 = pc 2 and 
T° l = cj l , respectively. Moreover, it is useful to define the GEM potentials $ and A 
in terms of p and j as hoo = 4$/c 2 and Jiqi = —2Ai/c 2 , respectively. Assuming that 
the source consists of a finite distribution of slowly moving matter with |v| << c, 
Ty ~ pviVj +p5ij, where p is the pressure, and (|1.3j) imply that hij = 0(c~ 4 ). 
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All terms of 0(c 4 ) will be neglected in this analysis. Under these conditions, the 
spacetime metric has the GEM form 

ds 2 = -c 2 (l - dt 2 - ^(A • dx)dt + (l + 2^ <%(&W. (1.4) 

In the Newtonian limit, $ reduces to the Newtonian gravitational potential, while 
A = 0(c _1 ). If the source distribution is confined around the origin of spatial 
coordinates, then far from the source 

, GM GJxx 

$ ~ , A~ =— , (1.5) 

r or 13 

where r = |x| and M and J are the mass and angular momentum of the source, 
respectively. Moreover, the transverse gauge condition reduces to 

;» +v (2 A ) =a (L6) 

We define the GEM fields via 

E = -V$-~ (\k\ , B = VxA, (1.7) 



cdt\2 

in direct analogy with electromagnetism. It follows from these definitions that the 
GEM fields have dimensions of acceleration and 

vxe =4IG b )- v G B ) =a (18> 

Furthermore, (|1.2j) implies that 

V-E = 4vrG /0 , Vx (^b) =i^E + — j. (1.9) 
V 2 / c at c 

The GEM field equations (|1.8|) and (jl.9p contain the continuity equation V • j + 
dp/dt = 0, as expected. 

For a complete GEM theory, we need an analogue of the Lorentz force law. The 
Lagrangian for the motion of a test particle of mass m, L = —mcds/dt, can be 
written to linear order in $ and A as 

(v 2 \ ^ / v 2 \ 2m 
1-^2 J + m l[ l + -^j^ ^7v-A, (1.10) 

where 7 is the Lorentz factor. The equation of motion, dp/dt = F, where p = 7mv 
is the kinetic momentum, takes a simple familiar form if dA/dt = and F is 
expressed to lowest order in v/c, <3? and A; then, 

F = -mE-2m-xB. (1.11) 

c 
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The canonical momentum of the particle is given in this case by p + (— 2m/c)A. 

Let us now discuss the gauge freedom of the GEM potentials. Under a coordinate 
transformation x^ — » x M = x^ 1 — e M , h^ u — > /i' = fo^,, + + e„ iAt to linear order in 
e^. Therefore, 

Ku = V + + e ",» - V^V (1-12) 

Under this transformation, the Riemann tensor remains invariant, but the connec- 
tion changes. We must restrict in such a way that those elements of the connection 
defining GEM fields also remain invariant. The new metric perturbation satisfies 
the transverse gauge condition as well, provided De^ = 0. Therefore, let eo = 0(c~ 3 ) 
and €i = 0(c -4 ); then, (|1.12p implies that 

$' = ®--jr^> ^U'^A + W, (1.13) 
cot 2 2 

where \I/ = c 2 e°/4 and = 0. Under the GEM gauge transformation (|1.13p . the 
GEM fields are invariant in close analogy with electrodynamics. 

It is important to discuss the stress-energy tensor for GEM in the context of 
our approximation scheme. For the sake of simplicity we set c = G = 1 in what 
follows, unless specified otherwise. The Landau-Lifshitz pseudotensor t^ u can be 
employed to determine the local stress-energy content of the GEM fields. The gen- 
eral gauge-dependent result is somewhat complicated and is given in [ID] : however, 
for a stationary configuration (i.e. d<fr/dt = and dA/dt = 0), one finds 

4vrGt o = ~\e 2 + H,i) A (i,fr (1-14) 
4irGt 0i = 2(E x B) is (1.15) 
AnGtij = [EiEj - ^E 2 ^ + (b^ + ^B^j . (1.16) 

There is some similarity between these and the corresponding relations in classical 
electrodynamics. In particular, the GEM Poynting vector is given by 

E x B. (1.17) 



2vrG 



For instance, gravitational energy circulates around a stationary source of mass m 
and angular momentum J = Jz with a flow velocity 

v, = ^sin^ (1.18) 

in the same sense as the rotation of the mass. Here we employ spherical polar coordi- 
nates and k = 4/7. The flow given by (|1.18p is divergence- free and the corresponding 
circulation is independent of the radial distance r and is given by 2nk{ J/M) sin 2 9. 
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1.3 Gravitational Larmor Theorem 

There is a certain degree of arbitrariness in the definitions of the GEM potentials and 
fields that leads to extraneous numerical factors as compared to electrodynamics. To 
develop GEM in a way that would provide the closest possible connection with the 
standard formulas of electrodynamics, we have adopted a convention that, among 
other things, provides a clear statement of the gravitational Larmor theorem. We 
assume that a test particle of inertial mass m has gravitoelectric charge qE = —m 
and gravitomagnetic charge qs = —2m. For a source that is a rotating body of mass 
M, the corresponding charges are positive, Qe = M and Qb = 2M, respectively, 
in order to preserve the attractive nature of gravity. Moreover, the ratio of the 
gravitomagnetic charge to the gravitoelectric charge is always 2, since linearized 
gravity is a spin-2 field. This is consistent with the fact that for a spin-1 field such 
as in Maxwell's theory, the corresponding ratio is always unity. 

The Larmor theorem originally established a basic local equivalence between 
magnetism and rotation [llj. In fact, the electromagnetic force on a test particle 
of mass m and charge q in the linear approximation is the same as the inertial 
force experienced by the free particle with respect to an accelerated system of ref- 
erence with translational acceleration sll = —qE^/m and frequency of rotation 

= 9_BB/(2mc). In electromagnetism qE = Qb = Q and for all particles with the 
same charge-to-mass ratio q/m, the electromagnetic field can be replaced by the 
same accelerated system. This circumstance takes on a universal character in the 
case of gravity since the gravitational charge-to-mass ratio is the same for all par- 
ticles according to the principle of equivalence of gravitational and inertial masses. 
The universality of the gravitational interaction thus leads to a geometric theory 
of gravitation, i.e. general relativity. An analogous approach to electrodynamics 
is impossible due to the fact that q/m for different particles can, for instance, be 
positive, negative or zero. 

To develop a gravitational analogue of Larmor's theorem, we recall from the 
previous section that in the linear approximation of general relativity, the exterior 
gravity of a rotating source can be described in terms of GEM fields. In a sufficiently 
small neighborhood of the exterior region, the GEM fields may be considered locally 
uniform. These fields may then be locally replaced by an accelerated system in 
Minkowski spacetime. To this end, let us imagine an accelerated observer following 
a worldline Xq(t). Here r is the observer's proper time and = dx^/dr and 
A^ = du^/dr are its velocity and acceleration vectors, respectively. Let 'VV) be the 
orthonormal tetrad frame of the observer such that = U>M and 

- 1 ^ = K^ {pV (1-19) 

where <j)ap{ T ) is the antisymmetric acceleration tensor of the observer. In analogy 
with the Faraday tensor, cft a p consists of an "electric" part 4>oi = cii and a "magnetic" 
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part 4>ij = eijkQ . Here a and are spacetime scalars that represent respectively 
the translational acceleration, a% = a^A^, and the rotational frequency of the local 
spatial frame with respect to the local nonrotating (i.e. Fermi- Walker transported) 
frame. Consider now a geodesic system of coordinates established along the 
worldline of the fiducial observer. At any event r along the worldline, the straight 
spacelike geodesic lines orthogonal to the worldline span a hyperplane that is Eu- 
clidean space. Let x^ be the coordinates of a point on this hyperplane; then, 

/ = xS + nyr), r = x°. (1.20) 

The Minkowski metric rj^dx^dx 1 " with respect to the new coordinates takes the 
form g^ u dX^dX u , where 

9 oo = -(l + a-X) 2 + (ft xX) 2 , (1.21) 
g 0i = (O X X);, gij = Sij. (1.22) 

These geodesic coordinates are admissible if goo < 0; a detailed discussion of the 
nature of the boundary of the admissible region is given in |12j . 

A comparison of the metric given by (jl.2ip and (jl.22p with (jl.4p reveals that 
they are Larmor equivalent at the linear order once 

* = -a x • X, A = ~{l L x X, (1.23) 

and we neglect spatial curvature. The corresponding GEM fields to lowest order are 
E = — V<3? = 3ll and B = V x A = — fi^, as expected from the traditional Larmor 
theorem with qE = —m and qs = — 2m. 

The gravitational Larmor theorem [13] is essentially Einstein's principle of equiv- 
alence formulated within the GEM framework. Einstein's heuristic principle of 
equivalence traditionally refers to the Einstein "elevator" and its translational ac- 
celeration in connection with the gravitoelectric field of the source. However, it 
follows from the gravitational Larmor theorem that a rotation of the elevator is 
generally necessary as well in order to take due account of the gravitomagnetic field 
of the source. 

In classical electrodynamics, a charged spinning test particle has a magnetic 
dipole moment n = qS/(2mc), where m, q and S are respectively the mass, charge 
and the spin of the particle. In an external magnetic field B, the test dipole has an 
interaction energy — fi ■ B and precesses due to a torque /ixB. In a similar way, a 
test gyroscope of spin S with q — > qs = —2m has a gravitomagnetic dipole moment 
fi g = — S/c and precesses in the exterior field of a rotating source of mass M and 
spin J with the frequency [U] 



«P = ^[3(J-r)f-J], 



(1.24) 
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where fip = B/c and the gravitomagnetic field B is given by the curl of the vector 
potential in (|1.5p that corresponds to a source with gravitomagnetic dipole moment 
J/c. The related interaction energy is S • ftp. A major aim of the GP-B is to 
measure (|1.24p for gyros in a polar orbit about the Earth. 

It follows from (jl.24p that 2ir/£lp is a characteristic timescale for the gravitomag- 
netic field. More generally, gravitomagnetic effects reveal an interesting temporal 
structure around a rotating mass; this can be further illustrated by the phenomena 
associated with the gravitomagnetic clock effect [15] and the gravitomagnetic time 
delay [H]. 

A more exact long-term post-Schwarzschild analysis of the orbital motion of 
an ideal test gyroscope in the field of a rotating source reveals that besides the 
gravitoelectric geodetic (i.e. de Sitter-Fokker) precession of the gyro axis there is 
a complex gravitomagnetic component involving precessional as well as nutational 
motions — the latter is known as relativistic nutation [17] ■ The net gravitomagnetic 
spin motion reduces in the post-Newtonian approximation to equation (11.241) . 

1.4 Spacetime Curvature Approach to GEM 

The main elements that underlie the gravitational Larmor theorem apply equally 
well but in a different context to an alternative treatment of GEM based on space- 
time curvature. Unlike the previous treatment, the new approach is not limited to 
perturbations of flat spacetime and can be employed in an arbitrary curved space- 
time. 

Consider a congruence of test observers following geodesies in a gravitational 
field. Choosing a reference observer in this congruence, we set up a Fermi coordinate 
system along its path. This amounts to constructing an inertial system of coordi- 
nates in the immediate neighborhood of the reference observer |18] . Let -^^(r) be 
the orthonormal tetrad of the reference observer. Here is the vector tangent 
to the worldline of the observer and r is the proper time along its path. The spatial 
frame i = 1>2,3, consists of unit vectors along ideal gyro directions that are 

parallel transported along the worldline. The Fermi frame is a geodesic reference 
system that is based on the nonrotating orthonormal tetrad ^r a y The metric of 
spacetime in Fermi coordinates = (T, X) is then given by 

g 00 = -1 - /»'„,(,; .V'.V + . . . , (1.25) 
50, = -^RojikXi X k + . . . , ( L26 ) 

9ij = Sij - ^R ik3l X k X l + . . . , (1.27) 

where R a ^^{T) is the projection of the Riemann curvature tensor on the orthonor- 
mal tetrad of the reference observer 

Ra/3f5 = Rnupa^a^lp^^^ls)- (I- 28 ) 
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Along the reference geodesic T = r, X = and = r]^ v by construction. The 
Fermi coordinates are admissible within a cylindrical spacetime region of radius ~ 1Z 
around the worldline of the reference observer. Here 1Z is the radius of curvature of 
spacetime. 

The metric in Fermi coordinates — within the limited region of admissibility — 
has the form of a perturbation about Minkowski spacetime; therefore, using the 
previous GEM approach and comparing (jl.25p and (|1.26p with (|1.4p . we define the 
new GEM potentials as 

$(T, X) = -i J R 0i0i (T)X i X^' + . . . , (1.29) 

A t {T, X) = ±Ro jik (T)XSx k + ..., (1.30) 

where the spatial curvature has been ignored. The GEM fields are defined in terms 
of the potentials as before and are given to lowest order as 



E i (T,X)=R 0iQj (T)X3 +..., (1.31 
1 
2 



Bi(T,X) = -^e ijk R jk ol (T)X l + .... (1.32) 



It is interesting to note that in this approach the gravitoelectric field is directly 
connected with the "electric" components of the curvature tensor i?o«Oj an d the 
gravitomagnetic field is directly connected with the "magnetic" components of the 
curvature tensor Roijk [H]- It is possible to combine (|1.3ip and (|1.32p in the GEM 
Faraday tensor 

F aj3 = —R a p 0i X l (1.33) 

to linear order in X, where Fq{ = —Ei and = eijkB k . Then, Maxwell's equations 
P[a/3,j] = an d F a ^ g = 4:irJ a are satisfied in this case to lowest order in |X|/7£ with 

4vr J a (T, 0) = -Ro a = -8vrG (r 0a - ^q 0a T^ (1.34) 

along the reference trajectory in Fermi coordinates. Here we have used the gravita- 
tional field equations as well as the symmetries of the Riemann curvature tensor. 

The new approach to GEM naturally contains the analogue of the Lorentz force 
law. The motion of free test particles in the congruence relative to the reference 
particle at the spatial origin of Fermi coordinates can be expressed as 

^ + RwtjXi + 2R ikl0 V k Xi + (2 J R 0fej0 ^V fc 

2 2 \ ( L35 ' ) 

+ -R lk] iV k V l + -RokjiV^V 1 ) Xi = 0, 

valid to linear order in the separation X. This geodesic deviation equation is a 
generalized Jacobi equation [20] in which the rate of geodesic separation (i.e. the 
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relative velocity of the test particle) V = dX./dT is in general arbitrary (|V| < 1 
at X = 0). To linear order in velocity, one can show that (|1.35p takes the Lorentz 
form 

<i 2 X 

= q E B + q B V x B, (1.36) 

where qE = —m and qs = —2m as before. 

The stress-energy tensor in the new approach can be constructed essentially from 
the Faraday tensor (|1 .33|) as in Maxwell's theory, i.e. 

GT a(i = (f^F^ - ^g^F^F^ , (1.37) 

where an extra factor of G has been introduced due to dimensional considerations. 
From equations (|1.33|) and (11.370 . we find 

T ^ = ^G ( Ra ^ R \ ~ {t^mR 15 ^ X l XK (1.38) 

This tensor as well as the Faraday tensor (jl.33p vanishes along the worldline of the 
reference observer in the Fermi system; indeed, this is an immediate consequence of 
the inertial character of this system along the reference trajectory and a realization 
of Einstein's principle of equivalence. Thus this treatment depends on our choice of 
a reference observer and the corresponding Fermi coordinate system. 

To obtain a coordinate-independent measure of the stress-energy content of the 
gravitational field, we invoke the notion that the physical measurement of such a 
quantity requires an averaging process [21]. Starting from an event (T, 0) on the 
reference worldline, we average the tensor given in ([1.381) over a small sphere of 
radius eL, where < e << 1 and L is an invariant length scale that is characteristic 
of the source of the gravitational field under consideration. For instance, L could be 
GM/c 2 or, in the absence of such a scale, the Planck length. The quadratic nature 
of (11.381) in the spatial coordinates implies that the averaging involves 

(X l X j ) = k(eL) 2 6 i:j , (1.39) 

where k = 1/5 or 1/3 depending on whether the averaging involves the volume or 
the surface of the sphere, respectively. In either case, the constant k can be absorbed 
in the definition of L. Thus 

ke 2 L 2 - 

l " Ta ^ = "4^ r ^' 0<7A f«) A (^) A (o) A fo)' ( L4 °) 
where T^ upa {x) is the Bel tensor given by 

W(x) = UrmBJ? + Rrf*zRJ P c ) ~ - A g^R aPpi R a ^ ■ (1.41) 
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This tensor bears a certain similarity with the Maxwell stress-energy tensor and 
on this basis was first defined by Bel [22] for Einstein spaces in 1958. The Bel 
superenergy tensor is symmetric and trace-free in its first pair of indices and only 
symmetric in the second pair of indices. In a Ricci-flat spacetime, the Riemann 
tensor reduces to the Weyl conformal tensor C^ upa and the Bel tensor reduces to 
the completely symmetric and trace-free Bel-Robinson tensor T^upa given by 

T^pa = -(C^Cj^ + C^Cjf) - —9^gpaC a p^C al3lS . (1-42) 

An invariant average GEM stress-energy tensor of the gravitational field can 
thus be defined up to a constant positive multiplicative factor by 



L 2 - 

T (a)(f3) = ^ T ^Pv\a) X tf) X (0) X (0)> ( L43 ) 

where is symmetric and traceless. In the Ricci-flat case, — ► TJ*^^, 

where is the gravitational stress-energy tensor. This designation refers to the 

fact that in a Ricci-flat spacetime, the spatial components of the curvature tensor 
in (ll.27h . which were essentially ignored in our GEM analysis, are indeed basically 
given by its electric components; therefore, T? a involves all of the components 
of the spacetime curvature. The stress-energy tensors and T* have properties 
reminiscent of Maxwell's electrodynamics and have been discussed in detail in [23] . 
For instance, in the exterior of a rotating mass the gravitational Poynting flux based 
on the new approach has a flow velocity given by (j!.18f) with k = 3. 

A significant generalization of the concept of superenergy tensors has been de- 
veloped in [23] . 

The gravitomagnetic contribution to the spacetime curvature due to a rotating 
source involves subtle cumulative effects |25j that can be measured in principle via 
relativistic gravity gradiometry [26]. 



1.5 Spin-Rotation-Gravity Coupling 

An issue of fundamental interest is whether intrinsic spin is affected by a gravito- 
magnetic field in basically the same way as the classical spin of an ideal gyroscope. 
This question is related to the inertia of intrinsic spin. The description of physical 
states in the quantum theory is based upon the irreducible unitary representations 
of the inhomogeneous Lorentz group, which are characterized by means of mass and 
spin. The inertial properties of mass in moving frames of reference are already well 
known: for instance, via Coriolis, centrifugal and other mechanical effects, as well 
as their quantum mechanical counterparts. The inertial properties of intrinsic spin 
involve the phenomena associated with the spin-rotation-gravity coupling. 

The coupling of intrinsic spin with rotation reveals the rotational inertia of in- 
trinsic spin. That is, as a particle moves, its intrinsic spin keeps its aspect with 
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respect to an inertial frame; therefore, the spin appears to rotate with respect to a 
rotating observer. To this motion of spin corresponds, according to quantum me- 
chanics, a Hamiltonian H = —jfl ■ S. The general formula for the transformation 
of energy turns out to be £ ' = ~f(£ — h£lA4), where f2 is the frequency of rotation 
of the observer and hM is the component of the total angular momentum along 
the axis of rotation; that is, M = 0, ±1, ±2, ... for a scalar or a vector particle, 
while A4 =F \ = 0, ±1, ±2, . . . for a Dirac particle. This formula relates the energy 
of a quantum system measured by a rotating observer £' to measurements per- 
formed in a global inertial frame and can be written in the JWKB approximation 
as £' = 7(5 — fl • J), where J = rxp + Sis the total angular momentum. Thus, 
£' = 7(5 — v • p) — 7$7 • S, so that in the absence of intrinsic spin we recover the 
classical expression for the energy of a particle as measured in the rotating frame 
with v = fi x r. The spin-rotation coupling therefore involves an energy shift given 
by the Hamiltonian H = —jfl ■ S [27j . 

Observational evidence for such an energy shift in the case of fermions has been 
provided in certain high-precision experiments by way of a small frequency offset due 
to the coupling between the nuclear spin of mercury and the rotation of the Earth 
[28\ [29] ; moreover, a direct approach using neutron or atom interferometry has 
been proposed in [30J. For photons, helicity-rotation coupling has been confirmed 
to rather high accuracy using rotating GPS receivers [31] : moreover, experimental 
evidence exists for such a coupling in the microwave and optical regimes in terms of 
the frequency shift of polarized radiation [30J. The modifications of Doppler effect 
and aberration due to the coupling of photon spin with the rotation of the source 
and/or receiver have been the subject of recent studies [32], [33] . 

Let us now turn to the coupling of spin with gravitomagnetic fields; the spin- 
gravity coupling is naturally related to the spin-rotation coupling by way of Ein- 
stein's principle of equivalence. That is, starting from the spin-rotation Hamiltonian, 
the transformation fi — ► — ftp leads, according to the gravitational Larmor theorem, 
to the spin-gravity Hamiltonian. 

It follows from these ideas that in Earth-based experiments, to every Hamilto- 
nian we must add the spin-rotation- gravity interaction Hamiltonian 5H = — fi^ • 
S + Sip ■ S, where fi® and dp refer to the rotation frequency of the Earth and 
the corresponding gravitomagnetic precession frequency, respectively. Thus in the 
approximation under consideration here a particle with intrinsic spin behaves es- 
sentially like an ideal gyroscope. The energy difference corresponding to a spin-1/2 
particle polarized vertically up and down relative to the surface of the Earth can be 
estimated from MIq = 10 -19 eV and hflp = 10 _29 eV. The measurement of the latter 
term is beyond present capabilities by several orders of magnitude. In this connec- 
tion, however, we note that near Jupiter Mlp = 10~ eV, and therefore it is likely 
that with further improvements in magnetometer design, the spin-gravitomagnetic 
coupling could become measurable in a satellite in orbit near the surface of Jupiter 
in the foreseeable future [34J. It is important to recognize that such a relativistic 
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quantum gravitational effect, like all other gravitational effects, is subject to the 
whole mass-energy content of the universe. It follows that our treatment has been 
based on certain cosmological assumptions regarding the distribution of angular 
momentum in the universe; specifically, we have assumed that on the largest scales 
there is no preferred sense of rotation. Moreover, in 5H the spin-gravity coupling 
term has a gradient. Therefore, there exists a gravitomagnetic Stern-Gerlach force 
— V(fip • S) on a spinning particle that is independent of the its mass and hence 
violates the universality of the gravitational acceleration. The weight of a body 
thus depends on its spin, but the effect is too small to be directly measurable in the 
foreseeable future. It is interesting to note that the Stern-Gerlach force has an exact 
analogue in the classical Mathisson-Papapetrou spin-curvature force. The results of 
this section are in agreement with the consequences of Dirac-type wave equations 
in the gravitational field of a rotating mass [35] . 
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